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We construct a two-dimensional electric circuit model for creeping discharge. Two types of dis- 
charge, surface corona and surface leader, are modeled by a two-step function of conductance. 
Branched patterns of surface leaders surrounded by the surface corona appear in numerical simu- 
lation. The fractal dimension of branched discharge patterns is calculated by changing voltage and 
capacitance. We find that surface leaders often grow stepwise in time, as is observed in lightning 
leaders of thunder. 



I. INTRODUCTION 
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I Various branching patterns such as dendritic crystals, river networks and blood vessel networks have been intensive!; 
studied from a view point of fractals. [11] Discharge patterns such as lightning are also typical branching patterns. 
In creeping discharge, a strong voltage is applied at a point electrode on the surface of dielectric materials. Discharge 
patterns appear around the point electrode at the interface between the dielectric materials and the surrounding gas 
I— i' or liquid. In a suitable range of voltages, densely branched patterns are observed, that are called Lichtenberg figures. 
C/^ . Fractal analysis was performed for the Lichtenberg figures obtained in experiments by several groups. 0,01 Kebbabi 
Qh ' and Bernounal found that the fractal dimension decreases with increasing thickness of dielectric materials. 
(— I , The DLA model is a simple stochastic model that can generate a fractal branching pattern. 5] The fractal dimension 
• ^ is 1.71 in the DLA model. However, there are a variety of discharge patterns, which have various fractal dimensions 
\ Df. In typical patterns in lightning discharge, the fractal dimension is rather small, e.g., Df ~ 1.2. For densely 
'— '. branched Lichtenberg figures, the fractal dimension is Df ^ 1.8. The DLA model was generalized for application to 
^ ' the various discharge patterns. Niemeyer et al. proposed the 77 model in which the growth probability is proportional 
, to E^, where E is the electric field. @ Wiesmann and Zeller proposed another model in which the growth probability is 
' proportional to E—Ec for E > Ec, where Ec is the threshold value of discharge. Kuoershtokh et al. proposed another 
Ps| ; stochastic model in which there is a transition from a weak discharge (streamer) to a strong discharge (leader) . [§] We 

■ proposed a deterministic electric circuit model of branching patterns in the barrier discharge where an alternating 
'""^ I voltage is applied between dielectric electrodes. @ 

■ Various complicated dynamics are observed in the growth process of discharge patterns. For creeping discharge, 
growth dynamics was studied experimentally using a high-speed camera. When the electric voltage V is increased 

. over its threshold value, surface corona appears locally around the central electrode. The corona discharge is a weak 
' discharge that appears locally around a sharp electrode. When V is beyond a second threshold, surface leaders 
^ , appear, which make a branching pattern. A surface leader shows a stronger discharge than the surface corona. The 
surface corona grows further around the tips of surface leaders. The length of surface leaders increases with y.fioj 
When a surface leader reaches the outer electrode, a flashover occurs, and a large current fiows, leading to dielectric 
breakdown. 

The growth of surface leaders is not always smooth in time, but intermittent in time.[ll| A similar intermittent 
growth of the leaders is well known in lightning discharge. [2] When a leader goes down from a thunder cloud to the 
ground, it moves in steps of about 30 m with a pause of about 40 ms between steps. This leader is called the stepped 
leader. When the leader reaches the ground, a flashover occurs and a strong flash called a return stroke appears, which 
we observe as lightning. The origin of stepped leaders is not completely understood. There are several qualitative 
theories of such stepped leaders. Bruce pointed out the importance of the transition from a weak glow discharge to 
a strong arc discharge. [l2j Kumar and Nagabhushana proposed a simulation model of stepped leaders based on a 
complicated electric breakdown process. ^3] 

In this paper, we propose a simple deterministic electric circuit model composed of resistors and capacitors for the 
creeping discharge and perform numerical simulation. We have assumed a two-step function for the conductance, 
which exhibits a transition from a weak discharge to a strong discharge. We show that densely branched patterns 
appear when the capacitance is relatively large. We flnd an intermittent growth of the discharge pattern, which is 
similar to the stepped leader. 



FIG. 1: (a) Schematic figure of tfie conductance a as a function of electric potential difference AV^. (b) Schematic figure of an 
electric circuit on a triangular lattice. 



II. ELECTRIC CIRCUIT MODEL 



Our electric circuit model is composed of capacitors and resistors whose conductance changes with the voltage. The 
conductance a of the resistance is assumed to change stepwise as shown in Fig. 1(a). The conductance a increases 
rapidly owing to the ionization, when voltage is increased beyond its threshold value and the discharge sets in. Taking 
another transition between the surface corona to the surface leader into account, we assume that there is another 
threshold and the conductance a increases further at the second threshold. We study a dynamical system for a and 
the electric potential V . The time evolution of a is assumed to obey 

= ai-a, for |Ay| < Ki, 

- CT2 - (T, for Vcl < \AV\ < Vc2, 

at 

= <T3-fT, for |Ay|>K2, (1) 

at 

where /S.V is the electric potential difference, Vd and Vc2 are the two threshold values, cri^2,3 are the stepped values 
of conductance, and r is a time constant. We further assume that the threshold values depend on ct as Vc = Vcq — aa, 
which makes the transitions discontinuous transitions accompanying hysteresis. Hysteresis is often observed in the 
ionization transition or a transition between different types of discharge states. 

The resistors are set on edges of a triangular lattice as shown in Fig. 1(b). The conductance of the resistor between 
two sites (i,j) and is denoted as crij,i',j>, and the electric potential at the (i,j) site is denoted as Vij. The 

conductance aij^i'j' obeys eq. (1) for AV ~ Vij — Vi',j'- That is, 

'^ ^^^'dt ^ (^i-'^^ for < Kio - 

r ^'^'i',' .3 0-2 - for Vcio - acri,j,t',j' < \ AVijA\j' \ < Vc2o - OL'^i..j,v .,y , 



dt 

da,. 



dt 



0-3 - a, for \AV^JJ'J'\ > Vc20 - aai,j^i',j'. (2) 



Each site is connected to the ground with a capacitor of capacitance C as shown in Fig. 1(b). The capacitors represent 
the dielectric material on which the creeping discharge occurs. The electric potential V of the central (1+6) sites is 
set to be Vq and the electric potential outside of a hexagonal area of side length L is set to be 0. The electric potential 
Vi,j at site obeys 

from the conservation law of the electric charge. The summation is taken over the nearest neighbor sites of 
(j,j). We show several numerical results for L = 75,T4io = 0.16, ^420 — 0.66, cti = 0.032, CT2 = 32, (J3 — 800, t = 0.1, 
and a — 0.02. C and Vq are changed as control parameters. However, similar results were also obtained for different 
parameter values. 
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FIG. 2: Discharge patterns at (a) Vb — 0.45, (b) 0.6, and (c) 0.7 for C = 0.035. Discharged edges whose electric potential 
difference is beyond the first threshold are plotted using dots. The solid lines denote discharged edges whose electric potential 
difference is beyond the second threshold. The inner and outer electrodes are drawn with dashed lines. 
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FIG. 3: Discharge patterns at (a) C = 0.024 and (b) 0.054 for Vo = 0.7. (c) Double logarithmic plot of R{t) vs. N{t) for 
C = 0.032 and 0.054. Two straight lines denote iV ~ 7?^ **^ and iV ~ R^ '"^^. (d) Fractal dimension D of discharged patterns for 
several C's. 



III. BRANCHING PATTERNS 



Figure 2 shows discharge patterns at three different values of Vg for C = 0.035. The initial conditions are aij^i'j' = 
0.001 + 0.0001 X r, where r is a random number between and 1, and Vij = 0. The initial conditions are slightly 
random, but there is no randomness in system parameters. The corona discharge appears for Vq > 0.23 in our system. 
For Vq < 0.22, all ai,j^i',j' 's approach cri, because the electric potential difference is below the first threshold at all sites. 
Figure 2(a) shows a discharge pattern at Vq = 0.45. The inner and outer hexagonal electrodes are drawn using dashed 
lines. (The inner electrode is hardly seen in these plots because it is too small.) The discharge pattern stops growing 
at a finite time. In this discharge pattern, <Tij,i\j' takes a value of <72 in the corona region. Corona discharge is a 
discharge state that appears only around a sharp electrode. For Vq > 0.50, the electric potential difference goes beyond 
the second threshold, a stronger discharge state appears, that is, a few surface leaders grow from the inner electrode. 
The surface corona appears around the surface leaders. Figure 2(b) shows a discharge pattern at Vq = 0.6. A random 
discharge pattern appears even in the uniform system. This is probably because chaotic dynamics or dynamical 
instability is involved in our deterministic system, similarly to the interface instability in crystal growth 14]. The 
discharge pattern reaches the outer electrode with y = for Vq > 0.61 in this configuration. Figure 2(c) is a discharge 
pattern at Vq — 0.7, which represents a flashover. 

The degree of ramification of branched patterns increases, as the capacitance C is increased. Figures 3(a) and 3(b) 
show discharge patterns at (a) C = 0.024 and (b) 0.054 for Vq — 0.7. For the quantitative analysis of the branched 
patterns, we have calculated the radius of gyration as 



R{t)^^T.(^l + y?)/N{t), (4) 

where the summation is taken over the discharged sites whose electric potential difference is beyond the first threshold, 
Xi and Hi are the x and y coordinates of the discharged site, and N is the total number of discharged sites at time 
t. Figure 3(c) shows double logarithmic plots of the relation of R{t) and N{t) for C = 0.032 and 0.054. The fractal 
dimension is evaluated from the slope of the relation of In. N{t) vs In. R{t) as J? ~ 1.38 for C = 0.024, and as D ^ 1.85 
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FIG. 4: (a) Discharge pattern under a uniform initial condition at C = 0.035 and Vb = 0.6. (b) Time evolution of the total 
number S{t) of edges included in the surface leaders. 
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FIG. 5: (a) Discharge pattern from a uniform initial condition at C = 0.035 and Vb ~ 0.6. The discharge occurs only on six 
straight lines. The second threshold is Vc20 = 0.64. (b) Time evolution of the total number S{t) of edges included in the surface 
leaders, (c) Time evolution of the electric potential V{i,j) on the central line j = L/2. (d) Time evolution of the conductance 
cri,j,i+i,j on the central line j — L/2. 

for C — 0.054. Figure 3(d) shows the fractal dimension as a function of C. As C is increased, the discharge pattern 
is more densely branched and the fractal dimension increases, although the discharge pattern stops finally and the 
fractal dimension cannot be precisely evaluated using a sufficiently large pattern. In the experiment by Kebbabi 
and Berounal, the fractal dimension becomes smaller as the thickness d of the dielectric material is increased. The 
capacitance of the dielectric material is inversely proportional to the thickness, i.e., C oc 1/d. Therefore, our numerical 
result is qualitatively consistent with their experimental result. 

IV. STEPPED LEADERS 

In the time evolution of the discharge pattern, the surface leaders and corona discharge interact with each other. 
That is, the surface leaders appear on some selected edges in the corona region, with a certain time delay after the 
appearance of corona discharge. The corona discharge further grows around the tips of the surface leaders. In the 
numerical simulations in the previous section, we have used a slightly random initial conditions for ct^j^^'j/'s. If 
the initial conditions are uniform, a symmetrical pattern appears. Figure 4(a) shows such a hexagonally symmetric 
discharge pattern at C = 0.035 and Vo — 0.6. The initial condition is (Jij^i\j> — 0.01 for all edges. Figure 4(b) shows 
the time evolution of the number S of edges where the electric potential difference is beyond the second threshold. 
The number S is proportional to the total length of the surface leaders. Figure 4(b) shows that the surface leaders 
grow stepwise in time. That is, the surface leaders stop and go intermittently. For example, the surface leaders grow 
rapidly at t ~ 0.45 from the inner electrode and stop a.tt^ 0.82. This behavior is similar to that of the stepped leaders 
in lightning discharge. Even under the previous random initial conditions, this type of stepped leader is observed for 
each branch, but the total number S{t) exhibits no such clear stepped behavior except for the first step from S* = 0, 
because the timing of the stop-and-go is different for each branch. 

To understand the behavior of the stepped leaders, we have simulated an even simpler model system where only 
the edges on the six straight lines can be discharged, as shown in Fig. 5(a). Figure 5(b) shows the time evolution of 
S{t) in this system for C = 0.035, Vq = 0.6, Vdo = 0.16, and Vc2o = 0.64. The intermittent behavior is observed 
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for 0.58 < 1^20 ^ 0.65, when other parameters are fixed. At 1420 = 0.66, the surface leaders do not grow in this 
system. For Vc2o < 0.57, the surface leaders move almost steadily. Figure 5 shows that the branching is not directly 
related to the origin of the stepped leaders, because no branching occurs in this simple system. Figure 5(c) shows a 
quasi-3D plot of the time evolution of V{i,j) on the central line j = L/2 for 20 < i < 80 between 1.2 < t < 2, and 
Fig. 5(d) show a quasi-3D plot of the time evolution of (Tij^i+ij on the central line j — L/2. The behavior of the 
stepped leaders is clearly observed in Fig. 5(d). The interaction between the corona discharge, surface leader, and 
electric potential is important for the intermittent motion. When Vc2o = Vcio, the difference between the two types 
of discharge disappears, and then we have observed no intermittent motion. These results are qualitatively consistent 
with previous observations by Shimazaki et al.|ll| and Bruce[l2|. 



V. SUMMARY AND DISCUSSION 



We have constructed a simple deterministic electric circuit model for creeping discharge. We have found that random 
branching patterns appear, even though the system parameters are uniform. We have found that more densely 
branched patterns appear as the capacitance C increases, which is qualitatively consistent with the experimental 
result. We have numerically found the behavior of stepped leaders, which seems to originate from the interaction with 
the corona discharge and surface leaders, although the detailed theoretical analysis of the stepped leaders is left for 
future study. 

A similar intermittent growth is observed in other systems. The stick-slip motion in the frictional oscillation is 
a typical example. The threshold by the maximum static friction plays an important role in the stick-slip motion. 
Wakita et al. found a periodic growth in a bacterial colony. T5l Some threshold values of bacterial density seem to play 
an important role in periodic growth. Nakanishi et al. found a periodic growth in the electrodeposition of Sn,[l^ and 
Sakaguchi et al. proposed a coupled map lattice model of the periodic growth. [l7j The transition from a facet growth 
to a diffusion-limited growth is induced by a certain threshold voltage in this electrodeposition system. We think that 
there is some analogy in these systems, in that there is a transition between different states via some threshold values. 
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